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We review the functional renormalization group (RG) approach to the BCS-BEC crossover for an 
ultracold gas of fermionic atoms. Formulated in terms of a scale-dependent effective action, the 
functional RG interpolates continuously between the atomic or molecular microphysics and the 
macroscopic physics on large length scales. We concentrate on the discussion of the phase diagram 
as a function of the scattering length and the temperature which is a paradigm example for the 
non-perturbative power of the functional RG. A systematic derivative expansion provides for both 
a description of the many-body physics and its expected universal features as well as an accurate 
account of the few-body physics and the associated BEC and BCS limits. 
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1. Introduction 

Many challenges in contemporary theoretical physics deal with strongly interacting quantum field 
theories or many-body systems. Progress often relies on the construction of exact or approximate 
solutions. In the absence of exact solutions, reliable and controlled approximation methods typically 
are the only source of information about the system and the underlying physical mechanisms. An 
approximation scheme may be considered as reliable and controlled if it is based on a systematic and 
consistent expansion scheme and shows convergence towards the exact result (which, however, is 
often not known). Textbook examples are, of course, provided by perturbative expansions or lattice 
discretizations both of which can be consistently evaluated to a given order or lattice refinement, 
systematically improved, and the convergence can at least be checked as a matter of practice. 

In this contribution, we would like to demonstrate that the functional RG can be used to develop 
systematic and consistent expansion schemes for strongly interacting systems. Most importantly, it 
can be applied in the spacetimc continuum and docs not require a perturbative ordering scheme. 
Nevertheless, it offers a variety of tools to verify qualitative and quantitative reliability and practical 
convergence. As a prime example of strongly interacting many-body systems, we take the BCS- 
BEC crossover as an illustration for the use of the functional RG. The concrete physical system 
that we have in mind is an ultracold atomic Fermi gas with two accessible hyperfine spin states 
near a Feshbach resonance, showing a smooth crossover between Bardeen-Coopcr-Schrieffcr (BCS) 
superfluidity and Bosc-Einstein condensation (BEC) of diatomic molecules (Eagles 1969; Leggett 
1980). 

By means of an external magnetic field B the phenomenon of a Feshbach resonance allows to 
arbitrarily regulate the effective interaction strength of the atoms, parameterized by the s-wave 
scattering length a. We briefly discuss the example of 6 Li (O'Hara et al. 2002), which is besides 40 K 
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realized in current experiments (Regal et al. 2004; Zwicrlein et al. 2004; Kinast et al. 2004; Bourdcl 
et al. 2004; Bartenstein et al. 2004; Partridge et al. 2005), see left panel of Fig. [TJ 

For magnetic fields larger than ~ f200G the scattering length a is small and negative, giv- 
ing rise to the many-body effect of Cooper-pairing and a BCS-type ground state below a critical 
temperature. The BCS ground state is superfluid described by a non-vanishing order parameter 
(f>o = (t/Jifa) bilinear in the fermion fields. An increase of the temperature leads to a second order 
phase transition to a normal fluid, (j>o = 0. Magnetic fields below B ~ 600G induce a small and 
positive scattering length a and the formation of a diatomic bound state, a dimer. The ground 
state is a BEC of repulsive dimers and again a phase transition from a superfluid, 0n > 0, to a 
normal fluid, cf>Q = 0, can be observed at a critical temperature. For magnetic fields in the regime 
700G < B < f 100G the modulus of the scattering length \a\ is large and diverges at the unitary 
point, B n = 834G where unitarity of the scattering matrix solely determines the two-body scatter- 
ing properties. At and near unitarity the fermions are in a strongly interacting regime. It connects 
the limits of BCS superfluidity and Bosc-Einstein condensation by a continuous crossover and also 
shows a superfluid ground state, with (j> > (Eagles 1969; Leggett 1980). 

A convenient parameterization of the crossover is given by the inverse concentration c _1 = 
(akp) 1 - Here the density of atoms n = kp/(3Tr 2 ) defines the Fermi momentum &f in natural 
units with h = ks = 2M = 1, where M is the mass of the atoms. The inverse concentration c _1 
varies from large negative values on the BCS side to large positive values on the BEC side with a 
zero-crossing at the unitary point, see right panel of Fig. [1] 

A description of the qualitative features of the BCS-BEC crossover has been achieved by Nozieres 
& Schmitt-Rink (1985) and Sa de Melo et al. (1993) within extended mean-field theories which 
account for the contribution of both fermionic and bosonic degrees of freedom. However, the quan- 
titatively precise understanding of BCS-BEC crossover physics requires non-perturbative methods. 
The experimental realization of molecule condensates and the subsequent crossover to a BCS-like 
state of weakly attractive fermions (Regal et al. 2004; Zwierlein et al. 2004; Kinast et al. 2004; 
Bourdel et al. 2004; Bartenstein et al. 2004; Partridge et al. 2005) pave the way to future exper- 
imental precision measurements and provide a testing ground for non-perturbative methods. An 
understanding of the crossover on a quantitative level at and near the resonance has been developed 
through numerical Quantum Monte-Carlo (QMC) methods (Carlson et al. 2003; Astrakharchik et 
al. 2004; Bulgac et al. 2006; Burovski et al. 2006; Akkincni et al. 2007). The complete phase di- 
agram has been accessed by functional field-theoretical techniques, as t-matrix approaches (Pieri 
& Strinati 2000; Perali et al. 2004), Dyson-Schwinger equations (Diehl & Wetterich 2006, 6), 2PI 




600 700 800 900 1000 1100 1200 Cooper pairs strongly interacting pairs diatomic molecules 

Magnetic field (G) 

Figure 1: Left: Schematic plot of a Feshbach resonance, see e.g. O'Hara (2002). Right: Sketch of the 
crossover physics. 
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methods (Haussmann et al. 2007), and RG flow equations (Birse et al. 2005; Dichl et al. 2007a, 
6; Gubbels & Stoof 2008; Bartosch et al. 2009; Krippa 2009). These pictures of the whole phase 
diagram (Nikolic & Sachdev 2007; Pieri & Strinati 2000, Perali et al. 2004; Diehl & Wetterich 2006, 
b; Birse et al. 2005; Diehl et al. 2007a, 6; Gubbels & Stoof 2008, Haussmann et al. 2007) do not yet 
reach a similar quantitative precision as anticipated for the QMC calculations. 

We intend to fill this gap and discuss the limit of broad Feshbach resonances for which all 
thermodynamic quantities can be expressed in terms of two dimensionless parameters, 



and the temperature in units of the Fermi temperature T/Tp where Tp = kp. In the broad resonance 
regime, macroscopic observables are universal (Diehl & Wetterich 2006, Nikolic & Sachdev 2007, 
Dichl et al. 2007a, Ho 2004), i.e. they are to a large extent independent of the concrete microscopic 
realization. Very similar to the notion of universality near 2nd-ordcr phase transitions, universality 
in the present context can be traced back to the existence of a fixed point in the RG flow which is 
approached provided the Feshbach (Yukawa) coupling is large enough (Diehl et al. 20076). 

Our review is based on the RG studies of (Diehl et al. 2007a, 6; Diehl et al. 2010; Floerchinger 
et al. 2009; Floerchinger et al. 2010). For RG studies using different expansion schemes, see (Birse 
et al. 2005; Gubbels k Stoof 2008; Bartosch et al. 2009; Krippa 2009). In the following, we first 
introduce the required techniques from RG flow equations for cold atoms, see Sees. -[b] Further, 
we include the quantitative effect of particle-hole fluctuations, Sec. QJ and systematically extend 
the truncation scheme, accounting for changes of the Fermi surface due to fluctuation effects, Sec. 
[HI Additionally, we also consider an atom-dimcr interaction term. 



Microscopically the BCS-BEC crossover can be described by an action including a two-component 
Grassmann field ip = (ipi,ip2), describing non-relativistic fermions in two hyperfine states and a 
complex scalar field 4> a s the bosonic degrees of freedom. In different regimes of the crossover (j) can 
be seen as a field describing molecules, Cooper pairs or simply an auxiliary field. Explicitly, the 
microscopic action at the ultraviolet scale A reads 



In thermal equilibrium, the system is described by the Matsubara formalism. The variable p, is 
the chemical potential. The parameter v/y(B) = v(B) + Sv(A) includes the detuning from the 
Feshbach resonance v(B) = pu(B — B ), with the magnetic moment of the boson field pm , and 
a renormalization counter term 5is(A) that has to be adjusted to match the conditions from the 
physical vacuum. The Yukawa coupling h\ is related to the width of the Feshbach resonance. Note 
that the bosonic field 4> appears quadratically in Eq. (|2 . 1 [) so the functional integral over qb can be 
carried out, corresponding to an inverse Hubbard-Stratonovich transformation, and our model is 
equivalent to a purely fermionic theory with an interaction term 



the concentration: 



c = akF , 



(1.1) 



2. Microscopic model and the functional RG 



S = 





^,e&(Pl + P2)lpl{p' 1 )^l(pi)lp2(P2)' l p2(p2) S(p! + p 2 ~ p\ ~ P 2 ) ■ 



(2.2) 



Article submitted to Royal Society 



4 



M. M. Scherer, S. Floerchinger and H. Gies 



Here, p = (po,p) and the microscopic interaction between the fermions is described by the tree-level 
expression with a classical inverse boson propagator in the denominator 

W«) = 5T-^ ^7 . ( 2 -3) 

-u>+\ -2fi + v A (B) 

where uj is the real-time frequency of the exchanged boson <p. It is related to the Matsubara frequency 
go via analytic continuation uj = —iqo. Further, q = pi + P2 is the center of mass momentum of the 
scattering fermions ipi and tp2 with momenta p\ and P2, respectively. The limit of broad Feshbach 
resonances corresponds to h\ — > oo. In this limit the microscopic interaction becomes pointlike, 
with strength —h A /u\. 

The functional RG can be formulated as a functional differential equation for an action functional 
for which the microscopic model serves as an initial value. Whereas the microscopic interactions 
are governed by S at the ultraviolet scale A, quantum and thermal fluctuations effectively modify 
the interactions at larger length scales which can be summarized in an effective action (e.g. for 
the IPI proper vertices) valid at a momentum scale k. In other words, includes the effects of 
fluctuations with momenta higher than k and governs the interactions with momenta near k. This 
effective average action or flowing action satisfies the Wetterich equation (Wetterich 1993), being 
an exact RG flow equation, 



d k r k m = isTr 



(rf[$] + i? fc ) 1 d k R k 



(2.4) 



Here, the STr operation involves an integration over momenta and a summation over internal 
indices with appropriate minus signs for fermions. The collective field $ summarizes all bosonic 

(2) 

and fcrmionic degrees of freedom, and 1^ [$] denotes the second functional derivative of T k 

(^H^^ra^w (2 - 5) 

The long-wavelength regulator R k specifies the details of the regularization scheme. Specific exam- 
ples will be discussed below. For reviews of the functional renormalization group see (Salmhofer & 
Honerkamp 2001; Berges et al. 2002; Gies 2006; Pawlowski 2007; Kopietz 2010). From the full ef- 
fective action in the long wavelength limit = rfc = o [<£>], all macroscopic properties of the system 
under consideration can be read off. 

Equation (|2.4j) is the technical starting point of our investigations. It is a functional differential 
equation which, upon expansion of this functional into a suitable basis translates to a system of 
infinitely many coupled differential equations for the expansion coefficients, i.e. generalized running 
couplings. Identifying suitable expansion schemes is not a formal but a physics problem: expan- 
sions should be based on building blocks that encode the relevant degrees of freedom of the system 
possibly at all scales. In the present context, this emphasizes the usefulness of composite bosonic 
fields which are expected to be the relevant long-range degrees of freedom at low temperatures. 
Reducing the full effective action to a treatable selection of generalized couplings defines a trunca- 
tion. Possible truncation schemes include vertex expansions, derivative expansions or other schemes 
to systematically classify all possible operators of a given system. The quantitative success of a 
given truncation scheme does not necessarily rely on the existence of a small expansion parameter 
like the interaction strength, but only requires that the operators neglected in a truncation do 



Article submitted to Royal Society 



Functional renormalization for the BCS-BEC crossover 



5 



not take a strong influence on the now of the operators included. In practice, a truncation can be 
tested in various ways, e.g., by verifying the practical convergence for increasing truncations or by 
studying regulator-scheme independencies for universal quantities. In the present context, also the 
comparison with well-known few-body physics turns out to provide a useful benchmark. 

3. Basic truncation 

(a) Derivative expansion 

Thermodynamics of a system can be obtained from its grand canonical partition function Z 
or the corresponding grand canonical potential VLq = —ThxZ. It is related to the effective action 
via r[$ eg ] = flc/T, where <E> e<? is obtained from the field equation ^r[$]|$ = $ e5 = 0. Let us first 
present a basic version of a truncation which already captures all the qualitative features of the 
BCS-BEC crossover: 

r*[$] = J ^(d T - A - fj,)tp + 0* {z^d T - ^) 4> + u(p, p) - h{4>*^2 + WMtVt- (3- 1 ) 

The effective potential U(p,p) is a function of p = (f>*(j) and p. This truncation can be motivated 
by a systematic derivative expansion and an analysis of the symmetries encoded in Ward identities 
(Dichl et al. 2007a; Diehl et al. 2010). It does not yet incorporate, for instance, the effects of 
particle-hole fluctuations and we will come back to this issue in Sect. |4] We define rcnormalizcd 

1/2 — - — 

fields <ft = 4>, p = A^p and renormalized couplings Z^ = Z^/A^, h = h/^A^ and express Eq. 
(|3.1[) in these quantities 

r*[S] = J - A - M + r [z^t - y) + U( P ,p) - h (<^iV>2 + <HZ4>i)\- (3-2) 

We expand the effective potential around the /c-dependent location of the minimum po(k) and the 
fc-independent value of the chemical potential pg corresponding to the physical particle number 
density n. We determine po(k) and po by the requirements (d p U)(po(k), p Q ) = for all k, and 
— (d fJl U)(po, po) = n at k = 0. More explicitly, we employ a simple expansion for U(p, p) of the form 

U(p,p) = U(p ,po) - n k (p - po) + (m 2 + a(p - p ))(p - po) + ^Hp ~ Po) 2 - (3-3) 

In the symmetric or normal gas phase, we have po = 0, while in the phase with spontaneous breaking 
of U(l) symmetry (superfluid phase), we have po > and m 2 = 0. The atom density n = —dU/dp 
corresponds to n k in the limit k — > 0. 

The running couplings in this truncation explicitly are m 2 (fc), A(fc), a(k), rife, Z^,(k) and h(k). 
In the phase with spontaneous symmetry breaking m 2 is traded for pg. In addition, we need the 
anomalous dimension r\ = — kd^hxA^. At the microscopic scale k = A the initial values of our 
couplings are determined from Eq. (|2.1[) . This gives m 2 (A) = V\{B) — 2p, po(A) = 0, A(A) = 0, 
Z^{K) = 1, h(A) = h\, a(A) = —2 and ua = 3ir 2 p6(p). Finally, our regularization scheme is 
specified by a regulator for space-like momenta which for the fermionic and bosonic field components 
reads 

R k ^ = (sign(/ - p)k 2 - {f - p)) 6 (k 2 - \f - p\) , R k j = A (k 2 - f/2) 6 (k 2 - f/2) , 
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respectively. For the fermions, it regularizes fluctuations around the Fermi surface, whereas bosonic 
fluctuations are suppressed for generic small momenta. This choice is optimized in the spirit of 
Litim (2000) and Pawlowski (2007). 

For our choice of the regulator and with the basic approximation scheme Eq. (|3.2[) the flow 
equation for the effective potential can be computed straightforwardly: 



\/2k 5 

kd k U = 7]P U' + -^r- 



3tt 2 Z, 

with the threshold functions 



k 



^((/i + k 2 )H(p + k 2 ) - (p - k 2 )H(p - k 2 ))s^, (3.4) 



s (0) - 

„(°) _ 



k 2 + U' 



k 2 +U'+2pU" 



, / k 2 + U'+2pU 77 \ ( 1 ■ AT 
+ V k?+U' )\2 +NE 



Vk 2 + U'y/k 2 +U'+2pU" 



y/k 



(3.5) 
(3.6) 



The threshold functions exhibit a temperature dependence via the Bose and Fermi functions 
•Wb/fIc] = (e £ / T =F I) -1 - From the effective potential flow, we derive the flow equations for the 
running couplings m 2 or po and A. For details we refer to (Dichl et at 2010). Further we need flow 
equations for Aa, and Z$ that are obtained by the projections 



oqo 



and d t A d 



«o=0 



2d t -^(P 4> ) 22 (0,q) 



q=0 



with the momentum-dependent part of the propagator 



5<fH(q)8<f> j (q') 



bi= \/2po-02— 



(3.7) 



(3.8) 



Here the boson field is expressed in a basis of real fields 4>(x) = -^(cf>i(x) + ifoix)). These flow 
equations arc derived by Diehl et al. (2010) and have a rather involved structure. Finally, we need 
the flow of the Yukawa coupling. In the symmetric regime the loop contribution vanishes and the 
flow is given by the anomalous dimension, 



dth — —rjh, or in dimensionless units dth 



(-l+ V )h 2 , 



(3.9) 



where h 2 = h 2 /k. In the regime of spontaneous symmetry breaking (po > 0) there is a loop 
contribution ~ h 3 \po from a mixed diagram involving both fermions and bosons. This contribution 
is quantitatively subleading which we have verified also numerically. For the basic approximation 
scheme, Eq. ()3.2|) . we therefore dropped this contribution. 



(b) Vacuum limit and contact to experiment 

The vacuum limit allows us to make contact with experiment. We find that for n = T = the 
crossover at finite density turns into a second-order phase transition in vacuum (Diehl & Wetterich 
2007, Nikolic & Sachdev 2007) as a function of the initial value to 2 (A). In order to see this, we 
consider the momentum-independent parts in both the fcrmion and the boson propagator, — p, 
(the "chemical potential" for the fermions in vacuum) and m(k = 0) , which act as gaps for the 
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propagation of fcrmions and bosons. We find the following constraints, separating two different 
branches of the physical vacuum (Diehl & Wetterich 2007), 

m 2 (0) > 0, ^ = atom phase (a -1 < 0), 

m 2 (0) =0, n < molecule phase (a" 1 > 0), (3.10) 
m 2 (0) =0, fj, = resonance (a -1 = 0). 

The initial values m 2 (A) and Ha can be connected to the two-particle scattering in vacuum close to 
a Fcshbach resonance. For this purpose, one follows the flow of m 2 (k) and h(k) in vacuum, e.g. on 
the BCS (a -1 < 0), i.e. \i = T = n = 0, and extracts the renormalized parameters m 2 = m 2 {k = 0), 
h = h(k = 0). They have to match the physical conditions formulated in Eq. (|3.10p . We obtain the 
two relations 

2/1 . . N h\ , h 2 (k = 0) h 2 (A) 
rh 2 (A) = m (B-B )-2fi+-±A, a = \ '— = K > 3.11 

07^ «7T TO^(fc = 0) »7T [lM(B — Bo) 

where /xm is the relative magnetic moment of the molecules. These relations fix the initial conditions 
of our model completely and similar reasoning confirms their validity on the BEC side. Now we can 
express the parameters m 2 (A) and h 2 (A) by the experimentally accessible quantities B — Bo and 
a. They remain valid also for non-vanishing density and temperature, as long as the UV cutoff A is 
much larger than T and p. 



(c) Many-body phase diagram 

Although our flow equations describe accurately the vacuum limit and can be used to determine 
interesting few-body parameters they are not restricted to that limit. In fact, for nonzero temper- 
ature and density, the flow deviates from its vacuum form at scales with k 2 < T or k 2 < Tp. The 
resulting system of ordinary coupled differential equations is then solved numerically for different 
chemical potentials (i and temperatures T. For temperatures sufficiently small compared to the 
Fermi temperature, T/Tp <C 1, we find that the effective potential U at the macroscopic scale 
k = develops a minimum at a nonzero field value po > 0, d p U(po) = 0. The system is then in the 
supcrfluid phase. For larger temperatures we find that the minimum is at po = and that the "mass 
parameter" m 2 is positive, m 2 — d p U(0) > 0. The critical temperature T c of this phase transition 
between the superfluid and the normal phase is then defined as the temperature where 

p =0, d p U{0) = at k = 0. (3.12) 

Throughout the whole crossover the transition po — > is continuous as a function of T demonstrat- 
ing that the phase transition is of second order. An analysis of the scaling of the correlation length 
confirms that the phase transition is governed by a Wilson-Fisher fixed point for the N = 2 univer- 
sality class throughout the crossover (Diehl et al. 2010). This reflects the fact that the symmetries 
are properly encoded also on the level of the truncated action. 

From the flow equations together with the initial conditions we can already recover all the 
qualitative features of the BCS-BEC crossover, e.g. compute the phase diagram for the phase 
transition to superfluidity. The result for this basic approximation is displayed in the right panel of 
Fig. [6] by the dot-dashed line. 
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(d) Fixed-point and universality 

In the vacuum limit and in the regime where k 2 ^> — fi the flow of the anomalous dimension reads 
n = h 2 /(6TT 2 k) (Diehl et al. 2010). Together with the dimensionless flow of the Yukawa coupling, 
Eq. (|3.9p . this reveals the existence of an IR attractive fixed point given by r\ = 1, h 2 = 6ir 2 . 
This fixed-point is approached rapidly if the initial value of h 2 (A)/A is large enough, i.e., in the 
broad- resonance limit. Then the memory of the microscopic value of h 2 (A) /A is lost at large length 
scales. Also all other parameters except for the mass term m 2 are attracted to IR fixed points, 
giving rise to universality. The fixed-point structure remains similar for non- vanishing density and 
temperature and these findings also apply in this regime and determine the critical physics of these 
non-relativistic quantum fields. For a given temperature, this fixed point has only one relevant 
direction which is related to the detuning of the resonance B — Bq. 

4. Particle-hole fluctuations 

(a) Gorkov's correction to BCS theory 

For small and negative scattering length c _1 < 0, |c| <C 1 (BCS side), the system can be treated 
by the perturbative BCS theory of superfluidity (Cooper 1956; Bardecn et al. 1957). However, 
there is a significant decrease of the critical temperature as compared to the original BCS result 
due to a screening effect of particlc-holc fluctuations in the medium (Gorkov & Mclik-Barkhudarov 
1961; Hciselberg et al. 2000). Here we will sketch the technique to include the effect of particle-hole 
fluctuations in our functional RG treatment as developed by Floerchinger et al. (2009). 

In an RG setting, the features of BCS theory can be described in a purely fermionic language with 
the fermion interaction vertex A^, as the only scale-dependent object. In general, the interaction 
vertex is momentum dependent, \^{p' l ,pi,p l 2 ,p2), and its flow has two contributions which are 
depicted in Fig. [2j including the external momentum labels. For fc — > 0, /i — !► 0, T — !► and n — > 
this coupling is related to the scattering length, a = -^\^(pi = 0). 

In the BCS approximation only the first diagram in Fig. [5J the particle-particle (pp) loop, 
is kept and the momentum dependence of the four fermion coupling is neglected, by replacing 
^■4>(p'iiP1tP2iP2) by the pointlikc coupling evaluated at zero momentum. For fi > 0, its effect 
increases as the temperature T is lowered. For small temperatures T < X^bcs the logarithmic 
divergence leads to the appearance of pairing, as X^, — > oo. In terms of the scattering length a, 
Fermi momentum hp and Fermi temperature Tf, the critical temperature is found to be 

T c ,bcs w0.61 T F e^/^M. (4.1) 
P{ 

Pi P2 

Figure 2: Running of the momentum-dependent vertex A^. Here dk indicates scale-derivatives with 
respect to the regulator in the propagators but does not act on the vertices. 
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At zero temperature the expression for the second diagram in Fig. [U the particle-hole (ph) loop, 
vanishes if it is evaluated for vanishing external momenta, as both poles of the frequency integration 
are always either in the upper or lower half of the complex plane. The dominant part of the 
scattering in a fermion gas occurs, however, for momenta on the Fermi surface rather than for 
zero momentum. For non-zero momenta of the external particles the particle-hole loop makes an 
important contribution. Setting the external frequencies to zero, we find that the inverse propagators 
in the particle-hole loop are 

Pl>(<l) = iqo + (q -Vi) 2 - V, and Pip{q) — iq + (<f — p'^) 2 — P- ( 4 -2) 

Depending on the value of the momenta pi and p%, there are now values of the loop momentum q 
for which the poles of the frequency integration arc in different half-planes so that there is a nonzero 
contribution even for T = 0. 

To include the effect of particlc-holc fluctuations one could take the full momentum-dependence 
of the vertex A^ into account. However, the resulting intcgro-diffcrcntial equations represent a 
substantial numerical challenge. As a simple and efficient approximation, one therefore restricts the 
flow to the running of a single coupling A^, by choosing an appropriate momentum projection. In 
the purely fermionic description, this flow equation has a simple structure and the solution for A^ 1 
can be written as 

(A v ,(fc = 0))" 1 = (A^fc = A))" 1 + pp-loop + ph-loop . (4.3) 

Since the ph-loop depends only weakly on the temperature, one can evaluate it at T — and add 
it to the initial value X 1 p(k = A) -1 . As T c depends exponentially on the "effective microscopic 

coupling" (A^ A ^ = (A^(ft = A) -1 + ph-loop), any shift in results in a multiplicative 

factor for T c . The numerical value of the ph-loop and therefore of the correction factor for T c /Tp 
depends on the precise projection description. The averaging prescription used by Gorkov & Mclik- 
Barkhudarov (1961) leads to 

Tc = (4 e )i/3 Tc ' BCS ~ ^ T <=>bcs (4.4) 
and similar for the gap A at zero temperature. 

(b) Scale- dependent bosonization 

In Sect. [2] we describe an effective four-fcrmion interaction by the exchange of a boson. In this 
picture the phase transition to the supcrfluid phase is indicated by the vanishing of the bosonic 
"mass term" m 2 = 0. Negative m 2 leads to the spontaneous breaking of U(l) symmetry, since 
the minimum of the effective potential occurs for a nonvanishing supcrfluid density po > 0. For 
m? > we can solve the field equation for the boson <j> as a functional of ip and insert the solution 
into the effective action. This leads to an effective four-fermion vertex describing the scattering 
^1(^1)^2(^2) -> ^1(^1)^2(^2) 

-h 2 

Aw,, e ff = — : T-rz ■>• (^-5) 

l(Pl + P2)0 + 2 (Pi + PV + m 

To investigate the breaking of U(l) symmetry and the onset of superfluidity, we first consider the 
flow of the bosonic propagator, which is mainly driven by the fermionic loop diagram. For the 
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effective four-fermion interaction this accounts for the particle-particle loop (see left panel, r.h.s. of 
Fig. [3]). In the BCS limit of a large microscopic m\ the running of m 2 for k — > reproduces the 
BCS result (Cooper 1956; Bardeen et al. 1957). 

The particle-hole fluctuations are not accounted for by the renormalization of the boson prop- 
agator. Indeed, we have neglected so far that a four-fermion interaction term A^, in the effective 
action is generated by the flow. This holds even if the microscopic pointlike interaction is absorbed 
by a Hubbard-Stratonovich transformation into an effective boson exchange such that A,/, (A) = 0. 
The strength of the total interaction between fermions 

-h 2 

*(pi +P2)o + § (ft +P2) 2 +m 2 

adds A^, to the piece generated by boson exchange. In the partially bosonized formulation, the flow 
of is generated by the box diagrams depicted in the right panel of Fig. [3] A direct connection 
to the particle-hole diagrams of Fig. [5] can be established on the BCS side and in the microscopic 
regime: There the boson gap m 2 is large. In this case, the effective fcrmion interaction in Eq. 
(|4.6p becomes momentum independent, diagrammatically corresponding to a contracted bosonic 
propagator. The box diagram in Fig. [3] is then equivalent to the particle-hole loop investigated in 
the last section with the pointlike approximation A^cff — > f° r the fcrmion interaction vertex. 

In contrast to the particle-particle fluctuations (leading to SSB for decreasing T), the particle- 
hole fluctuations lead only to quantitative corrections and depend only weakly on temperature. This 
can be checked explicitly in the pointlike approximation, and holds not only in the BCS regime 
where T/fx <C 1, but also for moderate T/[i as realized at the critical temperature in the unitary 
regime. We therefore evaluate the box diagrams in Fig. [5] for zero temperature. We emphasize that 
a temperature dependence, resulting from the couplings parameterizing the boson propagator, is 
implicitly taken into account. For the external momenta we use an averaging on the Fermi surface 
similar to the one of Gorkov & Mclik-Barkhudarov (1961). For details see Floerchinger et al. (2009). 

After these preliminaries, we can now incorporate the effect of particle-hole fluctuations in 
the RG flow. In principle one could simply take A,/, as an additional coupling into account. How- 
ever, it is much more elegant to use a scale-dependent Hubbard-Stratonovich transformation (Gies 
& Wetterich 2001; Pawlowski 2007; Floerchinger & Wctterich 2009) which absorbes A^ into the 
Yukawa- type interaction with the bosons at every scale k. By construction, there is then no self 
interaction between the fermionic quasiparticles. The general procedure of "partial bosonization" is 
discussed in detail in Floerchinger et al. (2009). A slightly modified scheme based on the exact flow 
equation derived in (Floerchinger & Wetterich 2009) has been used in (Floerchinger et al. 2010). In 
that formulation one finds for the renormalized coupling m 2 in the symmetric regime an additional 

Pi P'2 , 

P2 P2 



pi p 2 pi p-j 

Figure 3: Left panel: Flow of the boson propagator. Right panel: Box diagram for the flow of the 
four-fermion interaction. 
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term reflecting the absorption of into the fcrmionic interaction induced by boson exchange, 

d t m 2 = d t m 2 \ m + — 9 t A^,| H g. (4.7) 

Here <9im 2 | Hg and 3tA^,| Hg denote the flow equations when the Hubbard-Stratonovich transforma- 
tion is kept fixed. Since A,/, remains now zero during the flow, the effective four-fermion interaction 
A^, e ff is purely given by the boson exchange. However, the contribution of the particle-hole exchange 
diagrams is incorporated via the second term in Eq. (|4.7[) . The flow equation of all other couplings 
are the same as with fixed Hubbard-Stratonovich transformation. In the regime with spontaneous 
symmetry breaking we use 

dth = dth\ KS H — j^-dtX^ | Hg , 
dtPo = ^*Po| HS — 2-^9 t A^,| Hg , 

d t X = d t X\ m +2^d t X^\ m . (4.8) 

We emphasize that our non-perturbative flow equations go beyond the treatment by Gorkov & 
Mclik-Barkhudarov (1961) which includes the particle-hole diagrams only in a perturbative way. 
Furthermore, the inner bosonic lines h 2 / P<p(q) in the box diagrams include the center-of-mass mo- 
mentum dependence of the four-fermion vertex. This is neglected in Gorkov's pointlike treatment, 
and thus represents a further improvement of the classic calculation. Actually, this momentum de- 
pendence becomes substantial away from the BCS regime where the physics of the bosonic bound 
state sets in. The continuous description of dynamically transmuting degrees of freedom is a par- 
ticular strength of an RG description, as exemplified also in the context of QCD (Gics & Wetterich 
2004, Braun 2009). 



5. Running Fermion sector 

In this section we aim at a systematic extension of the truncation scheme and consider a running 
fermion sector. Similar paramcterizations of the fcrmionic self-energy have been studied in Gubbels 
& Stoof (2008), Bartosch et al. (2009) and Strack et al. (2008). Further, we include an atom-dimcr 
interaction term. This section is based on the work by Floerchingcr et al. (2010). 

(a) Completion of the truncation 

In addition to the running couplings that have occurred so far in Sees. El now we want to take 
into account fc-dependent parameters m 2 and Z$ , in order to parametrize fluctuation effects on the 
self-energy of the fcrmionic quasiparticles. The extension of the truncation explicitcly reads 

T k = ^ dr Jd 3 x^Z^(d T - A)^ + ml$H + (z^d T - 4> 

+u(p,fx) - M<f + H>li>l) + W*<^ j. (5.1) 

The additional inclusion of the atom-dimcr coupling A^ closes the truncation on the level of 
interaction terms quartic in the fields and describes three-body scattering (Diehl et al. 2007c). It 
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-8 -6 -4 -2 2 4 6 -2 -1 1 2 3 4 



Figure 4: Left panel: Flow of the inverse fermionic wave function renormalization l/Z^ at T = 
at three different points of the crossover: c _1 = —1 (dashed line), c _1 = (solid line), c _1 = 1 
(short dashed line). Right panel: Inverse fermionic wave function renormalization 1/Z^, at k = as 
a function of the crossover parameter c" 1 . 



leads to quantitative modifications for the many-body problem. In the regime of spontaneously 
broken symmetry (pq > 0) the atom dimer coupling leads to a modification of the Fermi surface, 
in addition to the gap y/ h 2 po. 

1 /2 — 1 /2 — 

We define the renormalized fields <f> = A^ <f>, p = A^p, ip = ip and study the flow of the 

rcnormalizcd couplings Z$ = Z^/A^, h = h/(A^ 2 Z^), X^ = X^ / (A^Z^), m| = rfi^/Z^. As 
before, we expand the effective potential in monomials of p, see Eq. (|3.3|) . We use again a purely 
space-like regulator which is adjusted to the running Fermi surface, 

Rk,* = Z^ [signup 2 - r%)k 2 - {f - r 2 F ))]0 (fc 2 - \f - r%\) , i? M = A$ [k 2 - f /2] 9 (k 2 - f /2) , 

where r 2 F = —to 2 — X^po- For the fermions it regularizes fluctuations around the running Fermi 
surface, while for the bosons fluctuations with small momenta are suppressed. 

As a first quantity we investigate the vacuum dimer-dimer scattering length cim expressed in 
units of the atom-atom scattering length a. On the BEC side of the resonance we can derive this 
quantity from the corresponding couplings by the equation 

— — = 2—^ — , A^.off = 8na = —7==, for p < and k = 0. (5.2) 

To explicitely compute the vacuum quantity ajvf/a, we choose a value for a on the far BEC side, 
where for broad resonances the identity a = (— /i) -1 / 2 holds. We evolve the flow of the couplings 
to the infrared and extract the value of A, completely fixing ajv//a. In this truncation including the 
Xfy^p vertex, we find om/b = 0.59, which is in very good agreement with the well-known result from 
a direct solution of the Schrodinger equation, aj\//a = 0.60 (Petrov et al 2004). The accuracy of this 
result is somewhat surprising since no momentum dependence of X^ has been taken into account. 
The latter has turned out to be important for the atom-dimer scattering Diehl et at (2007c). On 
the other hand, the general leading-order effect of fermionic momentum-dependencies is captured 
by the wave function renormalization Z^ the effect of which is largest in the strongly interacting 
regime, when |c _1 | < 1, cf. Fig. S) 
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Figure 5: Left panel: Effective Fermi radius rp /k F as a function of the crossover parameter c" 1 for 
vanishing temperature (solid line). We compare to the effective Fermi radius in an approximation 
without the contribution of the atom-dimer vertex A^y, (dotted line). Right panel: Positive branch 
of the dispersion relation u>(q) in units of Ef for c _1 = —1 (dashed), c _1 = (solid) and c _1 = 1 
(short dashed). 



(b) Fermi sphere and dispersion relation 

The dispersion relation can be computed from the determinant of the renormalized fermionic 
propagator 

q-1 _ ( -h<j> £ -u - (q 2 + m\ + X<wPo) 

^ \-oj + (g 2 + m| + X^po) h(j) e 

by the equation detG,^ 1 = 0. Here, we have evaluated G,^ 1 in the regime of spontaneously bro- 
ken symmetry and performed analytical continuation to real frequencies u>. We see that the dis- 
persion relation is affected by the running of the couplings Z^, and A^, which follows as 

lj = ±-\/A 2 + (q 1 — rp) 2 , where A = h^/po is the gap and rp = \J ~ m2 b ~ X^po is the effective 
radius of the Fermi sphere. On the far BCS side of the crossover the renormalization effects on 
rp are small and rp approaches its classical value rp — \fp — where kp = (S^n) 1 / 3 . Close 
to the resonance, the Fermi sphere gets smaller. It finally vanishes on the BEC side at a point 
with c" 1 « 0.6, see Fig. [5] Here, the fermions are gapped even for A — >• by a positive value of 




6. Results 

For the studies in section [SJ we have omitted the effect of particle- hole fluctuations for simplicity. 
In the following, however, all the results are given for the correspondingly extended truncation 
including particle-hole fluctuations. 



(a) Single-particle gap at T = 

As a first study including all the couplings introduced in this contribution we investigate the 
single particle gap at zero temperature. On the far BCS side it is possible to compare to the result 
by Gorkov & Mclik-Bakhudarov (1961) which is given by A/E F = (2/e) 7 / 3 e 7r /( 2c ). Our approach 
allows to extend to the strongly interacting regime and even to the BEC side of the crossover, see 
Fig- El At the unitary point, (akp)^ 1 = 0, we obtain A/Ep = 0.46. Further, we compare our result 
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A/E F 


Carlson et al. (2003) (QMC) 


0.43 


0.54 


Perali et al. (2004) (t-matrix approach) 


0.46 


0.53 


Haussmann et al. (2007) (2PI) 


0.36 


0.46 


Bartosch et al. (2009) (FRG, vertex exp.) 


0.32 


0.61 


Floerchinger et al. (2010) (FRG, derivative exp.) 


0.51 


0.46 



Table 1: Results for the single-particle gap and the chemical potential at T = and at the unitary 
point by various authors. 



for chemical potential in units of the Fermi energy at the unitary point fi/Ep = 0.51 to different 
(non-perturbative) methods in Tab. [TJ 



(6) Phase diagram 

Our results for the critical temperature for the phase transition to superfluidity throughout the 
crossover are shown in Fig. [5] We plot the critical temperature in units of the Fermi temperature 
T c /Tp as a function of the scattering length measured in units of the inverse Fermi momentum, 
i. e. the concentration c = akp. 

On the BCS side of the crossover, where c _1 <C — 1, the BCS approximation and the effect of 
particle-hole fluctuations yield a critical temperature (Gorkov & Melik-Barkhudarov 1961) 



Tp 71 



e T/(«M ~ 0.28e' r/(afci?) . (6.1) 



depicted by the short dashed line in the right panel of Fig. [51 Here, C « 0.577 is Euler's constant. 
On the BEC side for very large and positive c" 1 , our result approaches the critical temperature of 
a free Bose gas where the bosons have twice the mass of the fermions Mb = 2M. In our units the 
critical temperature is then 

^£ = « 0.218. (6.2) 

Tf (6tt 2 C(3/2)) 2/3 

In-between there is the unitarity regime, where the two-atom scattering length diverges (c _1 — > 0) 
and we deal with a system of strongly interacting fermions. 

Our best result including particle-hole fluctuations is given by the solid line. This may be 
compared with a functional rcnormalization flow investigation ignoring particle-hole fluctuations as 
discussed in section [3] (dot-dashed line) (Diehl et al. 2007a). For c — > 0_ the solid line of our result 
agrees with BCS theory including the correction by Gorkov & Melik-Barkhudarov (1961). Deviations 
from this perturbative regime appear only rather close to the regime of strong interactions c _1 — > 0. 

For c — > 0+ this value is approached in the form (Baym et al. 1999) 

T c - T C .BEC 1/3 O-M C 

— = na M n M =k , . (6.3) 

J c ,bec a (67T ) ' 

Here, um — n/2 is the density of molecules and a,M is the molecular scattering length. Using our 
result olm/o, = 0.59 obtained from solving the flow equations in vacuum, the coefficient determining 
the shift in T c compared to the free Bose gas yields k = 1.39, see also Diehl et al. (2007a). In (Arnold 
& Moore 2001; Kashurnikov et al. 2001), the result for an interacting BEC is determined as/t= 1.31 
(dashed curve on BEC side ofO, see also (Blaizot et al. 2006a, b) for a functional RG study. This is 
in reasonable agreement with our result. As further characteristic quantities we give the maximum 
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3 4 

Figure 6: Left panel: Gap in units of the Fermi energy A/Ep as a function of (akp) -1 (solid line). 
For comparison, we also plot the result found by Gorkov and Melik-Bakhudarov (dashed) and 
extrapolate this to the unitary point (akp)^ 1 — 0, where Aqmb/Ef = 0.49. Right panel: Critical 
temperature T c /Tp in units of the Fermi temperature as a function of the crossover parameter 
(akp) -1 . The solid line gives the result of our full computation refered to in the text. The dot- 
dashed line is obtained using the more basic truncation discussed in section[3l The dotted line shown 
for (akp) -1 < shows the result of the perturbative calculation by Gorkov & Melik-Barkhudarov 
(1961). The dashed line corresponds to an interacting BEC with the shift in T c according to Eq. 
(|6.3[) . We use here om/« = 0.60 and n = 1.31. The three red dots close to and at unitarity show 
the QMC results by Burovski et al. (2008), while the single purple dot gives the result of Akkineni 
et al. (2007). 



/j, c /Ep Tc/Tp 



Burovski et al. (2006) (QMC) 


0.49 


0.15 


Bulgac et al. (2006) (QMC) 


0.43 


< 0.15 


Akkineni et al. (2007) (QMC) 




0.245 


Previous FRG estimate by Floerchinger et al. (2009) 


0.68 


0.276 


Floerchinger et al. (2010) (FRG) 


0.55 


0.248 



Table 2: Results for T c /Tp and ii c /Tp at the unitary point by various authors. 



of the ratio (T c /Tp) maK rj 0.31 and the location of the maximum (akp)^ x ss 0.40. For c _1 > 0.5, 
the effect of the particle-hole fluctuations vanishes. This is expected, since the chemical potential 
is now negative y, < such that the Fermi surface disappears. 

In the unitary regime (c _1 0), the particle- hole fluctuations still have a quantitative effect. 
We can give an improved estimate for the critical temperature at the resonance (c _1 = 0) where we 
find T c /Tp = 0.248 and a chemical potential fi c /Tp = 0.55. A comparison to other methods and 
our previous work is given in Tab. [2] We observe reasonable agreement with QMC results for the 
chemical potential y c /Tp, however, our critical temperature T c /Tp is larger. 



7. Discussion and Outlook 

As illustrated with the example of the BCS-BEC crossover, the functional RG is capable of de- 
scribing strongly-interacting many-body systems in a consistent and controllable fashion. Once the 
relevant degrees of freedom are identified - possibly in a scale-dependent manner approximation 
schemes based on expansions of the effective action can be devised that facilitate systematically 
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improvable quantitative estimates of physical observables. For the BCS-BEC crossover, already a 
simple derivative expansion including fermionic and composite bosonic degrees of freedom exhibits 
all qualitative features of the phase diagram. 

The inclusion of particle-hole fluctuations and higher-orders of the derivative expansion improve 
our numerical results in the BCS as well as in the BEC limit of the crossover in agreement with 
well-known other field-theoretical methods. We obtain satisfactory quantitative precision on the 
BCS and BEC sides of the resonance. Remarkably, the functional RG allows for a description of 
both, many-body as well as few-body physics within the same formalism. For instance, our result 
for the molecular scattering length ratio Bm/o is in good agreement with the exact result (Petrov 
et al. 2004). This quantitative accuracy is remarkable, as we have started with a purely fermionic 
microscopic theory without propagating bosonic degrees of freedom or bosonic interactions. 

In the strongly interacting regime where the scattering length diverges, no exact analytical 
treatments arc available. Our results for the gap A/Ep and the chemical potential n/Ep at zero 
temperature are in reasonable agreement with Monte-Carlo simulations. This holds also for the 
ratio fic/Ep at the critical temperature. The critical temperature T c /Tp itself is found to be larger 
than the Monte-Carlo result. 

In future studies, our approximations might be improved mainly at two points. One is the 
frequency- and momentum dependence of the boson propagator. In the strongly interacting regime, 
this might be rather involved, developing structures beyond our current approximation. A more 
detailed resolution might lead to modifications in the contributions from bosonic fluctuations to 
various flow equations. Another point concerns structures in the fermion-fermion interaction that 
go beyond a diatom bound-state exchange process. Close to the unitary point, other contributions 
might arise, for example in form of a ferromagnetic channel. While further quantitative modifications 
in the unitarity regime are conceivable, the present approximation already allows for a coherent 
description of the BCS-BEC crossover for all values of the scattering length, temperature and density 
by one simple method and microscopic model. This includes the critical behavior of a second order 
phase transition as well as the vacuum, BEC and BCS limits. 

In order to improve the comparison between the QMC simulations and the functional RG 
results in the strongly interacting regime, the Wetterich equation can also be evaluated in a finite 
volume. This may shed light on possible finite size effects in the QMC simulations, and can help 
to quantitatively compare finite-volume studies with infinite-volume calculations inherent to most 
analytical works. 

The authors are grateful to J. Braun, S. Diehl, J. M. Pawlowski, C. Wetterich for collaboration on the 
subject reviewed here. This work has been supported by the DFG research unit FOR 723. H.G. acknowledges 
support by the DFG under contract Gi 328/5-1 (Heisenberg program). S.F. acknowledges support by the 
Helmholtz Alliance HA216/EMMI. 
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